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Abstract 

In the framework of the theory of open systems based on completely posi- 
tive quantum dynamical semigroups, we give a description of the continuous- 
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ly-^ , variable quantum entanglement and quantum discord for a system consisting 

(Nl | of two noninteracting modes embedded in a thermal environment. Entangle- 

ment and discord are used to quantify the quantum correlations of the system. 
For all values of the temperature of the thermal reservoir, an initial separable 
Gaussian state remains separable for all times. We study the time evolution 
of logarithmic negativity, which characterizes the degree of entanglement, and 
show that in the case of an entangled initial Gaussian state, entanglement sup- 
pression (entanglement sudden death) takes place, for non-zero temperatures 
of the environment. Only for a zero temperature of the thermal bath the ini- 
tial entangled state remains entangled for finite times. We analyze the time 
evolution of the Gaussian quantum discord, which is a measure of all quan- 
tum correlations in the bipartite state, including entanglement, and show that 
quantum discord decays asymptotically in time under the effect of the thermal 
bath. This is contrast with the sudden death of entanglement. Before the sup- 
pression of the entanglement, the qualitative evolution of quantum discord is 
very similar to that of the entanglement. We describe also the time evolution of 
the degree of classical correlations and of quantum mutual information, which 
measures the total correlations of the quantum system. 

1 Introduction 

The study of quantum correlations is a key issue in quantum information theory 
[1] and quantum entanglement represents the indispensable physical resource for the 
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description and performance of quantum information processing tasks, like quan- 
tum teleportation, cryptography, superdense coding and quantum computation [2]. 
However, entanglement does not describe all the non-classical properties of quantum 
correlations. For instance, recent theoretical and experimental results indicate that 
some non-entangled mixed states can improve performance in some quantum com- 
puting tasks [3]. Zurek pfl [5] defined the quantum discord as a measure of quantum 
correlations which includes entanglement of bipartite systems and it can also exist 
in separable states. Very recently, quantum discord was related with the problem 
of irreversibility of entanglement [6] and an operational interpretation was given to 
quantum discord in terms on consumption of entanglement in an extended quantum 
state merging protocol |8] . The total amount of correlations contained in a quan- 
tum state is given by the quantum mutual information which is equal to the sum of 
the quantum discord and classical correlations [9]. 

In recent years there is an increasing interest in using nonclassical entangled 
states of continuous variable systems in applications of quantum information process- 
ing, communication and computation [10J. A full characterization of the nonclassical 
properties of such states exists, at present, only for the class of Gaussian states. In 
this special case there exist necessary and sufficient criteria of entanglement [TTJ [12] 
and quantitative entanglement measures [131 [H]. 111 quantum information theory of 
continuous variable systems, Gaussian states, in particular two-mode Gaussian states, 
play a key role since they can be easily created and controlled experimentally. 

Implementation of quantum communication and computation encounters the dif- 
ficulty that any realistic quantum system cannot be isolated and it always has to inter- 
act with its environment. Quantum coherence and entanglement of quantum systems 
are inevitably influenced during their interaction with the external environment. As 
a result of the irreversible and uncontrollable phenomenon of quantum decoherence, 
the purity and entanglement of quantum states are in most cases degraded. Practi- 
cally, compared with the discrete variable entangled states, the continuous variable 
entangled states may be more efficient because they are less affected by decoherence. 

Due to the unavoidable interaction with the environment, any pure quantum 
state evolves into a mixed state and to describe realistically continuous variable 
quantum information processes it is necessary to take decoherence and dissipation 
into consideration. Decoherence and dynamics of quantum entanglement in continu- 
ous variable open systems have been intensively studied in the last years [I5j [161 El 
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In this work we study, in the framework of the theory of open systems based 
on completely positive quantum dynamical semigroups, the dynamics of the continu- 
ous variable quantum entanglement and quantum discord of a subsystem consisting 
of two uncoupled modes (harmonic oscillators) interacting with a common thermal 
environment. We are interested in discussing the correlation effect of the environ- 
ment, therefore we assume that the two modes are independent, i.e. they do not 
interact directly. The initial state of the subsystem is taken of Gaussian form and 
the evolution under the quantum dynamical semigroup assures the preservation in 
time of the Gaussian form of the state. We have studied previously [321 [33] the evo- 
lution of the entanglement of two identical harmonic oscillators interacting with a 
general environment, characterized by general diffusion and dissipation coefficients. 
We obtained that, depending on the values of these coefficients, the state keeps for 
all times its initial type: separable or entangled. In other cases, entanglement gener- 
ation, entanglement sudden death or a periodic collapse and revival of entanglement 
take place. In [M] it was analyzed the non-Markovian short-time-scale evolution of 
entanglement and quantum discord for initial two-mode squeezed thermal vacuum 
states of a system consisting of two identical and non-interacting harmonic oscillators 
coupled to either two independent bosonic baths or to a common bosonic bath. In the 
independent reservoirs case, it was observed the detrimental effect of the environment 
for these quantities. In the common reservoir case, for initial uncorrelated states, it 
was found that only quantum discord can be created via interaction with the bath, 
while entanglement remains absent. 

The paper is organized as follows. In Sec. 2 we write the Markovian master 
equation in the Heisenberg representation for two uncoupled harmonic oscillators 
interacting with a general environment and give the general solution of the evolution 
equation for the covariance matrix, i.e. we derive the variances and covariances of 
coordinates and momenta corresponding to a generic two-mode Gaussian state. By 
using the Peres-Simon necessary and sufficient condition for separability of two-mode 
Gaussian states [TlJ [35] , we investigate in Sec. 3 the dynamics of entanglement for the 
considered subsystem. For all values of the temperature of the thermal reservoir, an 
initial separable Gaussian state remains separable for all times. We analyze the time 
evolution of the logarithmic negativity, which characterizes the degree of entanglement 
of the quantum state, and show that in the case of an entangled initial Gaussian 
state, entanglement suppression (entanglement sudden death) takes place, for non- 
zero temperatures of the environment. Only for a zero temperature of the thermal 
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bath the initial entangled state remains entangled for all finite times, but in the 
limit of infinite time it evolves asymptotically to an equilibrium state which is always 
separable. We analyze the time evolution of the Gaussian quantum discord, which is 
a measure of all quantum correlations in the bipartite state, including entanglement, 
and show that quantum discord decays asymptotically in time under the effect of 
the thermal bath. This is contrast with the sudden death of entanglement. Before 
the suppression of the entanglement, the qualitative evolution of quantum discord is 
very similar to that of the entanglement. We describe also the time evolution of the 
degree of classical correlations and of quantum mutual information, which measures 
the total amount of correlations of the quantum system. A summary is given in Sec. 
4. 

2 Equations of motion of two independent modes 
interacting with the environment 

We study the dynamics of the subsystem composed of two non-interacting modes in 
weak interaction with a thermal environment. In the axiomatic formalism based on 
completely positive quantum dynamical semigroups, the irreversible time evolution 
of an open system is described in the Heisenberg representation by the following 
quantum Markovian Kossakowski-Lindblad master equation for an operator A (f 
denotes Hermitian conjugation) [36 | [371 138] : 

f = {[H,A] + IUVJ[A,V 3 ] + [Vj , A]Vj). (1) 

Here, H denotes the Hamiltonian of the open system and the operators Vj, Vj , defined 
on the Hilbert space of H, represent the interaction of the open system with the 
environment. 

We are interested in the set of Gaussian states, therefore we introduce such 
quantum dynamical semigroups that preserve this set during time evolution of the 
system and in this case our model represents a Gaussian noise channel. Consequently 
H is taken a polynomial of second degree in the coordinates x, y and momenta p x ,p y 
of the two quantum oscillators and Vj,Vj are taken polynomials of first degree in 
these canonical observables. Then in the linear space spanned by the coordinates and 
momenta there exist only four linearly independent operators ^=1,2,3,4 [39J : 

Vj = a xj p x + dyjPy + b xj x + b yj y, (2) 
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where a X j, a y j } b X j, b y j are complex coefficients. The Hamiltonian of the two uncoupled 
non-resonant harmonic oscillators of identical mass m and frequencies ui and LO2 is 



H 



m 



.2J2 



2„.2\ 



(3) 



The fact that the evolution is given by a dynamical semigroup implies the pos- 
itivity of the matrix formed by the scalar products of the four vectors sl x , h x , & y , h y , 

X j,b X j,a y j,b y j, respectively. We take this matrix 
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It follows that the principal minors of this matrix are positive or zero. From the 
Cauchy-Schwarz inequality the following relations hold for the coefficients defined in 
Eq. OH): 
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We introduce the following 4x4 bimodal covariance matrix: 
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The problem of solving the master equation for the operators in Heisenberg repre- 
sentation can be transformed into a problem of solving first-order in time, coupled 
linear differential equations for the covariance matrix elements. Namely, from Eq. 
(JT]) we obtain the following system of equations for the quantum correlations of the 
canonical observables (T denotes the transposed matrix) [39|: 

da(t) 



dt 



Ya(t)+a(t)Y L + 2D, 



(7) 
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The time-dependent solution of Eq. ([7]) is given by [39] 

a(t) = M(t)[a(0) - a(oo)]M T (t) + a(oo), (10) 

where the matrix M(t) = exp(Yt) has to fulfill the condition lim^oo M(t) = 0. In 
order that this limit exists, Y must only have eigenvalues with negative real parts. 
The values at infinity are obtained from the equation 

Ya(oo) + a(oo)Y T = -2D. (11) 



3 Dynamics of continuous variable entanglement 
and discord 

To describe the dynamics of quantum correlations, we use two types of measures: 
logarithmic negativity for entanglement, and quantum discord. 

3.1 Time evolution of entanglement and logarithmic nega- 
tivity 

A well-known sufficient condition for inseparability is the so-called Peres-Horodecki 
criterion [35, 40J, which is based on the observation that the non-completely positive 
nature of the partial transposition operation of the density matrix for a bipartite sys- 
tem (transposition with respect to degrees of freedom of one subsystem only) may turn 
an inseparable state into a nonphysical state. The signature of this non-physicality, 
and thus of quantum entanglement, is the appearance of a negative eigenvalue in the 
eigenspectrum of the partially transposed density matrix of a bipartite system. The 
characterization of the separability of continuous variable states using second-order 
moments of quadrature operators was given in Refs. [HI [12]. For Gaussian states, 
whose statistical properties are fully characterized by just second-order moments, 
this criterion was proven to be necessary and sufficient: A Gaussian continuous vari- 
able state is separable if and only if the partial transpose of its density matrix is 
non- negative (positive partial transpose (PPT) criterion). 
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The two-mode Gaussian state is entirely specified by its covariance matrix (J6]) 
which is a real, symmetric and positive matrix with the following block structure: 



= A « > ( 12 ) 




where A, 5 and C are 2x2 Hermitian matrices. A and S denote the symmetric 
covariance matrices for the individual reduced one-mode states, while the matrix C 
contains the cross-correlations between modes. When these correlations have non- 
zero values, then the states with det C > are separable states, but for det C < it 
may be possible that the states are entangled. 

The 4x4 covariance matrix ([12"]) (where all first moments can be set to zero 
by means of local unitary operations which do not affect the entanglement) contains 
four local symplectic invariants in form of the determinants of the block matrices 
A,B,C and covariance matrix a. Based on the above invariants, Simon [TT] derived 
the following PPT criterion for bipartite Gaussian continuous variable states: the 
necessary and sufficient condition for separability is S(t) > 0, where 

S(t) = detAdetB + (^- |detC|) 2 
— Tr [AJC JB JC T J] - -(det A + det B) (13) 
and J is the 2x2 symplectic matrix 

We suppose that the asymptotic state of the considered open system is a Gibbs 
state corresponding to two independent quantum harmonic oscillators in thermal equi- 
librium at temperature T. Then the quantum diffusion coefficients have the following 
form 1 381: 



1 1 



mui 2 2kT' 

mu 2 D yy = = ^ coth (15) 

mw 2 2 2kl 

D X p x — D yPy — D xy — Dp xPy — D xpy — Dyp x — 0. 

The elements of the covariance matrix can be calculated from Eqs. f llOp . fill]) - 
Solving for the time evolution of the covariance matrix elements, we can obtain the 
entanglement dynamics by using the Simon criterion. 
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For Gaussian states, the measures of entanglement of bipartite systems are based 
on some invariants constructed from the elements of the covariance matrix [TBI 1201 |2T] . 
In order to quantify the degree of entanglement of the infinite-dimensional bipartite 
system states of the two oscillators it is suitable to use the logarithmic negativity. 
For a Gaussian density operator, the logarithmic negativity is completely defined by 
the symplectic spectrum of the partial transpose of the covariance matrix. It is given 
by En = — log 2 2z>_, where z/_ is the smallest of the two symplectic eigenvalues of the 
partial transpose a of the 2- mode covariance matrix a [20] : 



2v- 



,~,2 



A =|= y A 2 — 4 det a 



(16) 



and A is the symplectic invariant (seralian) , given by A = det A + det B — 2 det C. 
In our model, the logarithmic negativity is calculated as 



E N (t) = --log 2 [4g{a(t))}, 



(17) 



where 



g(a(t)) = -(det A + det B) - det C 

2 \ 1/2 

- det a{t) 



- (det A + det B) - det C 



It determines the strength of entanglement for E^{t) > 0, and if E^(t) < 0, then the 
state is separable. 

In the following, we analyze the dependence of the Simon function S(t) and of 
the logarithmic negativity E^it) on time t and temperature T of the thermal bath, 
with the diffusion coefficients given by Eqs. ( TT5l) . We consider two types of the initial 
Gaussian states: 1) separable and 2) entangled. 

1) We consider a separable initial Gaussian state, with the two modes initially 
prepared in their single-mode squeezed states (unimodal squeezed state) and with its 
initial covariance matrix taken of the form 



*.(0) = g 



/ cosh r 
sinhr 


V o 



sinhr \ 

cosh r 

cosh r sinh r 

sinh r cosh r / 



(19) 



where r denotes the squeezing parameter. In this case S(t) becomes strictly positive 
after the initial moment of time (5(0) = 0), so that the initial separable state remains 
separable for all values of the temperature T and for all times. 
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2) The evolution of an entangled initial state is illustrated in Figure 1, where 
we represent the dependence of the logarithmic negativity E^(t) on time t and tem- 
perature T for an entangled initial Gaussian state, taken of the form of a two-mode 
vacuum squeezed state, with the initial covariance matrix given by 



*.(0) = \ 



I cosh r sinh r \ 

cosh r — sinh r 

sinhr coshr 

V — sinh r cosh r j 



(20) 



We observe that for a non-zero temperature T, at certain finite moment of time, 
which depends on T, E^if) becomes zero and therefore the state becomes separable. 
This is the so-called phenomenon of entanglement sudden death. It is in contrast 
to the quantum decoherence, during which the loss of quantum coherence is usually 
gradual [29l W\\ . For T = 0, E^it) remains strictly positive for finite times and tends 
asymptotically to for t — > oo. Therefore, only for zero temperature of the thermal 
bath the initial entangled state remains entangled for finite times and this state tends 
asymptotically to a separable one for infinitely large time. One can also show that the 
dissipation favorizes the phenomenon of entanglement sudden death - with increasing 
the dissipation parameter A, the entanglement suppression happens earlier. 

The dynamics of entanglement of the two oscillators depends strongly on the 
initial states and the coefficients describing the interaction of the system with the 
thermal environment (dissipation constant and temperature). As expected, the loga- 
rithmic negativity has a behaviour similar to that one of the Simon function in what 
concerns the characteristics of the state of being separable or entangled [32 | |33 | FHl W2\ . 

3.2 Asymptotic entanglement 

On general grounds, one expects that the effects of decoherence is dominant in the 
long-time regime, so that no quantum correlations (entanglement) is expected to be 
left at infinity. Indeed, using the diffusion coefficients given by Eqs. ffl5|) . we obtain 
from Eq. ( TlTi) the following elements of the asymptotic matrices A{oo) and B(oo) : 

muia xx {oo) = = -coth— — , a xpx {oo) = 0, 

mu 2 a yy {oo) = a ' p y p y^ = I co th a yPy (oo) = (21) 

and of the entanglement matrix C(oo) : 

0-^(00) = 0-^(00) = a yPx (oo) = a PxPy (oo) = 0. (22) 
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Figure 1: Logarithmic negativity versus time t and temperature T for an entangled 
initial vacuum squeezed state with squeezing parameter r = 4 and A = 0.1, uj\ — 002 — 
1. We take m = U — k — 1. 

Then the Simon expression (|T3l) takes the following form in the limit of large times: 

and, correspondingly, the equilibrium asymptotic state is always separable in the case 
of two non-interacting harmonic oscillators immersed in a common thermal reservoir. 

In Refs. [321 EH1 HH E2} H3J HI] we described the dependence of the logarithmic 
negativity E]y(t) on time and mixed diffusion coefficient for two harmonic oscillators 
interacting with a general environment. In the present case of a thermal bath, the 
asymptotic logarithmic negativity is given by (for uj\ < 002) 

^(oo) = -log 2 coth^. (24) 

It depends only on temperature, and does not depend on the initial Gaussian state. 
En (00) < for T 7^ and E^{oo) = for T = 0, and this confirms the previous 
statement that the asymptotic state is always separable. 

3.3 Gaussian quantum discord 

The separability of quantum states has often been described as a property synony- 
mous with the classicality. However, recent studies have shown that separable states, 
usually considered as being classically correlated, might also contain quantum correla- 
tions. Quantum discord was introduced (HE] as a measure of all quantum correlations 



10 



in a bipartite state, including - but not restricted to - entanglement. Quantum dis- 
cord has been denned as the difference between two quantum analogues of classically 
equivalent expression of the mutual information, which is a measure of total correla- 
tions in a quantum state. For pure entangled states quantum discord coincides with 
the entropy of entanglement. Quantum discord can be different from zero also for 
some mixed separable state and therefore the correlations in such separable states 
with positive discord are an indicator of quantumness. States with zero discord rep- 
resent essentially a classical probability distribution embedded in a quantum system. 

For an arbitrary bipartite state pu, the total correlations are expressed by quan- 
tum mutual information [15] 

I<Jto)= E S( Pi )-S(p l2 ), (25) 

i=l,2 

where pi represents the reduced density matrix of subsystem i and S(p) = — Tr(p In p) 
is the von Neumann entropy. Henderson and Vedral [9J proposed a measure of bipar- 
tite classical correlations C{p\ 2 ) based on a complete set of local projectors {n|} on 
the subsystem 2: the classical correlation in the bipartite quantum state p\ 2 can be 
given by 

C(p 12 ) = S( Pl ) -mf m {S(p ll2 )}, (26) 

where S(pi\ 2 ) = 12kP k S(Pi) is the conditional entropy of subsystem 1 and inf {S(pi\ 2 )} 
represents the minimal value of the entropy with respect to a complete set of local 
measurements {n|}. Here, p k is the measurement probability for the kth local pro- 
jector and p\ denotes the reduced state of subsystem 1 after the local measurements. 
Then the quantum discord is defined by 

D(p 12 ) = I(p 12 )-C(p 12 ). (27) 

Originally the quantum discord was defined and evaluated mainly for finite di- 
mensional systems. Very recently [4*6l W7\ the notion of discord has been extended to 
the domain of continuous variable systems, in particular to the analysis of bipartite 
systems described by two-mode Gaussian states. Closed formulas have been derived 
for bipartite thermal squeezed states jl6] and for all two-mode Gaussian states [47J. 

The Gaussian quantum discord of a general two-mode Gaussian state p\ 2 can 
be defined as the quantum discord where the conditional entropy is restricted to 
generalized Gaussian positive operator valued measurements (POVM) on the mode 
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2 and in terms of symplectic invariants it is given by (the symmetry between the two 
modes 1 and 2 is broken) |47] 

D = f(y/0) - f{v-) - f(v+) + f(y/e), (28) 

where 

v x + 1 x + 1 x-1 x-1 
f W = log — — log — ^— , (29) 



2 7 2 + (/3 - - a) + 2| 7 | v / 7 2 + (/3-l)(5-a) 
if (5 - a/3) 2 <(P + l) 7 2 (a + 5) 

(30) 

a/3 - 7 2 + 5 - + (5 - a/3) 2 - 2 7 2 (5 + a/3) 

2^ ' 

otherwise, 



a = 4detA, (3 = 4detB, 7 = 4detC, 5 = 16deta, (31) 

and z/ T are the symplectic eigenvalues of the state, given by 

2z/ 2 = A =f VA 2 - 4 det a, (32) 

where A = det A + det B + 2 det C. Notice that Gaussian quantum discord only 
depends on | detC|, i.e., entangled (det C < 0) and separable states are treated on 
equal footing. 

The evolution of the Gaussian quantum discord D is illustrated in Figure 2, 
where we represent the dependence of D on time t and temperature T for an entangled 
initial Gaussian state, taken of the form of a two-mode vacuum squeezed state (12"U|) . 
for such values of the parameters which satisfy for all times the first condition in 
formula fl30|) . The Gaussian discord has nonzero values for all finite times and this 
fact certifies the existence of nonclassical correlations in two-mode Gaussian states 
- either separable or entangled. Gaussian discord asymptotically decreases in time, 
compared to the case of the logarithmic negativity, which has an evolution leading to 
a sudden suppression of the entanglement. For entangled initial states the Gaussian 
discord remains strictly positive in time and in the limit of infinite time it tends 
asymptotically to zero, corresponding to the thermal product (separable) state, with 
no correlation at all. One can easily show that for a separable initial Gaussian state 
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Figure 2: Gaussian quantum discord D versus time t and temperature T for an 
entangled initial vacuum squeezed state with squeezing parameter r = 4 and A = 
0.1, u)\ — 0J2 — 1- We take m — % — k — 1. 



with covariance matrix (|T9|) the quantum discord is zero and it keeps this values 
during the whole time evolution of the state. 

From Figures 1 and 2 we notice that, in concordance with the general properties 
of the Gaussian quantum discord [47], the states can be either separable or entangled 
for D < 1 and all the states above the threshold D = 1 are entangled. We also notice 
that the decay of quantum discord is stronger when the temperature T is increasing. 

It should be remarked that the decay of quantum discord is very similar to that 
of the entanglement before the time of the sudden death of entanglement. In the 
vicinity of a zero logarithmic negativity (E^ = 0), the nonzero values of the discord 
can quantify the nonclassical correlations for separable mixed states and one considers 
that this fact could make possible some tasks in quantum computation [48J. 

The measure of classical correlations for a general two-mode Gaussian state p\2 
can also be calculated and it is given by [47] 



while the expression of the quantum mutual information, which measures the total 
correlations, is given by 



In Figures 3 and 4 we illustrate the evolution of the degree of classical correlations 
C and, respectively, of the quantum mutual information /, as functions of time t 




(33) 




(34) 
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Figure 3: Degree of classical correlations C versus time t and temperature T for 
an entangled initial vacuum squeezed state with squeezing parameter r = 4 and 
A = 0.1, oj\ — 002 — 1. We take m — h — k — 1. 

and temperature T for an entangled initial Gaussian state, taken of the form of a 
two- mode vacuum squeezed state (1201) . These two quantities manifest a qualitative 
behaviour similar to that one of the Gaussian discord: they have nonzero values for 
all finite times and in the limit of infinite time they tend asymptotically to zero, 
corresponding to the thermal product (separable) state, with no correlation at all. 
One can also see that the classical correlations and quantum mutual information 
decrease with increasing the temperature of the thermal bath. 

4 Summary 

In the framework of the theory of open systems based on completely positive quantum 
dynamical semigroups, we investigated the Markovian dynamics of quantum correla- 
tions for a subsystem composed of two noninteracting modes embedded in a thermal 
bath. We have presented and discussed the influence of the environment on the dy- 
namics of quantum entanglement and quantum discord for different initial states. We 
have described the time evolution of the logarithmic negativity, which characterizes 
the degree of entanglement of the quantum state, in terms of the covariance matrix for 
Gaussian input states, for the case when the asymptotic state of the considered open 
system is a Gibbs state corresponding to two independent quantum harmonic oscil- 
lators in thermal equilibrium. The dynamics of the quantum entanglement strongly 
depends on the initial states and the parameters characterizing the environment (dis- 
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Figure 4: Quantum mutual information / versus time t and temperature T for an 
entangled initial vacuum squeezed state with squeezing parameter r = 4 and A = 
0.1, u)\ — 0J2 — 1- We take m — % — k — 1. 

sipation coefficient and temperature). For all values of the temperature of the thermal 
reservoir, an initial separable Gaussian state remains separable for all times. In the 
case of an entangled initial Gaussian state, entanglement suppression (entanglement 
sudden death) takes place for non-zero temperatures of the environment. Only for a 
zero temperature of the thermal bath the initial entangled state remains entangled 
for finite times, but in the limit of infinite time it evolves asymptotically to an equilib- 
rium state which is always separable. The time when the entanglement is suppressed, 
decreases with increasing the temperature and dissipation. 

We described also the time evolution of the Gaussian quantum discord, which is 
a measure of all quantum correlations in the bipartite state, including entanglement. 
The values of quantum discord decrease asymptotically in time. This phenomenon 
is different from the sudden death of entanglement. The time evolution of quantum 
discord is very similar to that of the entanglement before the sudden suppression of 
the entanglement. After the sudden death of the entanglement, the nonzero values of 
quantum discord manifest the existence of quantum correlation for separable mixed 
states. Quantum discord is decreasing with increasing the temperature. One considers 
that the robustness of quantum discord could favorize the realization of scalable 
quantum computing in contrast to the fragility of the entanglement [48J . 

Presently there is a large debate relative to the physical interpretation existing 
behind the fascinating phenomena of quantum decoherence and existence of quantum 
correlations - quantum entanglement and quantum discord. Due to the increased 
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interest manifested towards the continuous variables approach [10J H9] to quantum 
information theory, the present results, in particular the existence of quantum discord 
and the possibility of maintaining a bipartite entanglement in a thermal environment 
for long times, might be useful in controlling entanglement and discord in open sys- 
tems and also for applications in the field of quantum information processing and 
communication. 
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